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ABSTRACT 

Aims. We study the properties and nature of extremely red galaxies (ERO, R — K >5.6) found behind two lensing clusters and compare 
them with other known galaxy populations. 

Methods. New HST/ACS observations, Spitzer IRAC and MIPS, and Chandra/ACIS observations of the two lensing clusters Abell 
1835 and AC1 14 contemplate our earlier optical and near-IR observations (Richard et al. 2006) and have been used to study extremely 
red objects (EROs) in these deep fields. 

Results. We have found 6 and 9 EROs in Abell 1835 and AC1 14. Several (7) of these objects are undetected up to the / and/or zsso 
band, and are hence "optical" drop-out sources. The photometric redshifts of most of our sources (80%) are z ~ 0.7-1.5. According 
to simple colour-colour diagrams the majority of our objects would be classified as hosting old stellar populations ("ellipticals"). 
However, there are clear signs of dusty starbursts for several among them. These objects correspond to the most extreme ones in R - K 
colour. We estimate a surface density of (0.97±0.31) arcmin 2 for EROs with (R-K >5.6) at K< 20.5. Among our 15 EROs 6 (40 %) 
also classify as distant red galaxies (DRGs). 11 of 13 EROs (85 %) with available IRAC photometry also fulfil the selection criteria 
for IRAC selected EROs (IEROs) of Yan et al. (2004). SED modelling shows that ~ 36 % of the IEROs in our sample are luminous 
or ultra-luminous infrared galaxies ((U)LIRG). Some very red DRGs are found to be very dusty starbursts, even (U)LIRGs, as also 
supported by their mid-IR photometry. No indication for AGNs is found, although faint activity cannot be excluded for all objects. 
From mid-IR and X-ray data 5 objects are clearly classified as starbursts. The derived properties are quite similar to those of DRGs 
and IEROs, except for 5 extreme objects in terms of colours, for which a very high extinction {Ay ^ 3) is found. 

Key words. Galaxies - high-redshift - evolution- starburst- Cosmology- early Universe- Infrared: galaxies 



1. Introduction 



Since their discovery in the late 1980ies ( Elston et al.|1988 1989 ), extremely red objects (EROs) have attracted serious attention 



These first detections were initially presumed to be high-redshift (z > 6) galaxies in a star-forming phase ( |Elston etal 
colour follow-up observations later identified these objects as luminous galaxies at z 



1988) . Multi- 



ar population 



0.8, dominated by an old stell; 

( |Elston etal.|1989) . The detection of two bright (K>18.4) extended objects (HR10 & HR14) with (I-K) colours near 6.5 by Hu and 
Ridgway (1994) highlighted the difficulty in classifying these galaxies. When first discovered, HR10 and HR14 were interpreted 
as being ellipticals at z ~ 2.4. Subsequent spectroscopic and morphological observations indicated that HR10 is not a quiescent 
elliptical galaxy, but rather a bright interacting galaxy at z = 1.44 (Graham & Dey 1996; |Stern et al.|2006[ ). 

In general there are two main scenarios which would produce a red enough spectral energy distribution to satisfy the established 
colour criteria for EROs (e.g. R-K >5-7, I-K >4-6), either due to an old passively evolved population or by extreme dust 
reddening as found in star bursts ( |Cowie et al.|1994| |Cimatti et al.|1999| |Daddi et al.|2002| |Geo rgakakis et al.|2006[) in a redshift 



range of 1 ^ z 55 2. A number of review articles discuss various aspects related to this topic (e.g. |McCarthy|20 04 ; Ferguson et al. 
12000] ). 

Especially the abundance of massive old ellipticals poses a strong test for the two competing scenarios of elliptical galaxy formation: 



early assembly (z/ >2-3) , e.g. by monolithic collapse, and passive luminosity evolution thereafter (PLE models) (Tinsley & Gunn 
1976 Pozzetti et al. 1996), or hierarchical merging of smaller sized objects ( White & Rees|1978 [Kauffmann et al. 1993| Somerville 



Base d on observations collected at the Very Large Telescope (Antu/UTl), Eur opean Southern Observatory, Paranal , Chile (ESO Programs 
69.A-0508, 70.A-0355, 73.A-0471), the NASA/ESA Hubble Space Telescope obtained at the Space Telescope Science Institute which is operated 



by AURA under NASA contract NAS5-26555, the Spitzer Space Telescope, which is operated by the Jet Propulsion Laboratory, California Institute 
of Technology under NASA contract 1407, and the Chandra satellite. 
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et al.||200"T) ). Observational evidence has been found for both scenarios: several surveys detected a deficit of ellipticals at z >1 



supporting the hierarchical merging models ( Roche et al.|2003 Kitzbichler & Wfiite|2 006), while others are consistent with PLE 
Pm et al.|2002l|Cimatfi et al.|2002[[Somerville et al.|2004) . 

However, in recent years the hierarchical merging scenario in a ACDM universe has been established as the favoured model. 
Nevertheless, the vast number of different renditions leaves room for dramatically different predictions regarding critical parameters 
like the number density of massive galaxies at specific times (Fontana et al.|2004 Treu et aLp 005 and references therein). 
The picture is complicated by the results of numerous morphological studies on EROs, which assigns a large fraction of EROs to 
disk galaxies at somewhat lower redshifts ( Yan & Thompson 2003 ; Gilb ank et al.|2003j Moustakas et al.|2004|>. In addition, a small 
fraction of EROs could also be act ive galactic nuclei (AGNs), as shown by deep XMM and Chandra data ( Alexander et al.|2 002 
|Roche et al.|2003l|Brusa et al.|2005| ). 



With the increasing number of large scale surveys like UKIDSS (Simps on et al.|20 06) and others it became clear that other means 
than spectroscopy is needed in order to classify to which of the two major galaxy populations the large number of EROs belong. 
These could either be combinations of RIJHK (Pozzetti & Mannucci 2000; Bergstrom & Wiklind 2004) colours, or the use of near 
and mid-infrared bands (Wilson et al. 2004). 

Independently of which exact colour criteria has been used, all EROs have at least one mutual property, their faintness at optical 
wavelengths causing limitations to the accuracy of photometric redshift estimates and other parameters derived from SED features. 
In order to increase the apparent brightness of EROs we use the natural magnification effect provided by massive galaxy clusters. 



This method has been applied successfully for the investigation of other faint sources, like Lyman break galaxies (|Pettini e t al. 2000 
ISwinbank et al.|2007) , faint SCUBA sources ( |Smail et al.|1998"l|Ivison et al.|2001| > and EROs ( [Smith et al.|2"002"l|Takata et al.|2003] > 



In the present study we proceed to a systematic selection of EROs in the fields of the two lensing clusters Abell 1835 and 
AC 114, based on observations obtained by Richard et al. (2006) and new ACS/HST, Spitzer, and Chandra observations. These are 
used to discuss their empirical properties, their nature and to derive physical parameters like photometric redshift, extinction, star 
formation rates, and stellar population properties. Several of these objects were found earlier in our //-band selected sample of 



optical drop-out objects constructed for the search of high redshift galaxies (see Richard et al. 2006| l. The ERO subsample from that 
paper is analysed in detail in Schaerer et al. (2007 1. 

Throughout this paper we adopted the following cosmology: Q m = 0.3, H o =70 km s 'Mpc -1 in a flat universe. All magnitudes 
are given in the Vega system if not stated otherwise. 



2. Observations and data reduction 



The observations described here are part of multi-colour observations on two galaxy clusters, AC 114 and Abell 1835, which have 
well known lensing properties. An extensive description of the initial observations (optical and near-infrared data) and the available 
data can be found in Richard et al. (2006). Exposure time, limiting magnitudes and more characteristics can be found in Tab. [T] 
The near-infrared data {SZ, J, H and K s ) were obtained with the Infrared Spectrometer and Array Camera (ISAAC) located on the 
VLT-UT1 (FOV 2.5 acrmin x 2.5 arcmin with a pixel size of 0.148 arcsec). The optical data for Abell 1835 (VRI) are archive data 
from the CFHT12k camera at CFHT (|Czoske et al.|2003|l, those for AC1 14 (UBVRI) were taken from|Campusano et al.|([200Ti. 



Table 1. Characteristics of the observations. Filters, effective wavelength, exposure time, depth (lcr), and AB corrections (Cab) 
correspond to Mab = '«vega + Cab are listed. a If only one filter name is given, than the same setup was used for both fields. h 
For small fractions of the final Rjq2 image, the total exposure time is only 400 seconds. For objects without R-band detection we 
adjusted the detection limits accordingly (see Tabj2]>. * The 24 //m integration time is as much as 3600 s within a ~30"wide strip 
crossing the cluster center because of the way multiple data sets were taken. 
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2.1. ACS data 

New z-band (F850LP, denoted z%so hereafter) imaging was obtained with the ACS camera onboard HST in November 2004 (AC1 14) 
and July 2005 (Abell 1835). The total observing time for AC114 and Abell 1835 were 9184 and 9110 seconds respectively. While 
the AC1 14 field is centred at the central cluster galaxy, Abell 1835 was observed off centre in order to avoid bright sources north of 
the cluster. For the reduction of both data sets we used the IRAF/Pyraf package multidrizzle on pos t-calibrated datapj 
For the calculation of the lcr detection limit, we applied the same method as for all other bands (see R ichard et al.|2006 i: measuring 



the corresponding standard deviation of the flux in randomly distributed circular apertures of 1.5 arcsec diameter (approx. 3000). 
None of the aperture positions lies within 3 arcsec to its closest neighbour or sources which have at least a 3cr detection. 

2.2. IRAC & MIPS data 



The 3.6, 4.5, 5.8, and 8.0 pm images were obtained using the Infrared Array Camera (IRAC; [Fazio et aL|2004 i while the 24 pm 



images were obtained using the Multi-band Imaging Photometer for Spitzer (MIPS; Rieke et al. 2004), both on board the Spitzer 
Space Telescope ( |Werner et al.|2004| ). The instruments, data, and reduction procedures are described in Egami et al. (2006 ). 



2.3. Chandra 

Both Abell 1835 and AC1 14 have been observed previously by Chandra. AC1 14 has been observed once for a total exposure of 75 
ksec (OBSID 1562). A comprehensive analysis of the cluster X-ray properties based on this dataset has been published previously 
in |De FilTp pis et alT] ( |2004| l. In the case of Abell 1835, two short archival observations are available (OBSIDs 495 and 496) as 
well as a deep 200 ksec GO observation (OBSIDs 6880, 6881, and 7370) obtained as part of this program. All available data were 
reprocessed using CIAO 3.2 and the latest calibration files available in CALDB 3.0. Standard screening was applied to all event 
files to remove bad grades, bad pixels, and background flares. After standard cleaning, the resulting net exposures were 73 and 224 
ksec for AC1 14 and Abell 1835, respectively. 

Bright X-ray point sources in the fields of AC114 and Abell 1835 were identified using the CIAO tool wavdetect. No bright 
X-ray point sources were detected within 5 arcsec of the ERO source positions. To determine flux upper limits, source spectra were 
extracted in a 2 arcsec radius aperture around each ERO source position. This aperture captures virtually all of the Chandra PSF 
over the range of relevant off-axis angles. The local background for each source was determined using an annular aperture from 
2-4 arcsec surrounding the source aperture. Count-weighted detector response (RMFs) and effective area (ARFs) files were created 
for each extraction region using the CIAO tools mkacisrmf and mkwarf, including the temporal, spectral, and spatial dependence 
of the ACIS filter contaminant. For Abel 1835, source spectra and matching background, RMF, and ARF files were produced for 
each OBSID separately and then fit jointly during the spectral analysis. This analysis is discussed in more detail in Section 3.3. All 



spectral analysis was done using the ISIS Houck & Denicola (2000) spectral fitting package and the XSPEC model library. 



In the case of Abell 1835, spectra from all 5 observations (2 archival datasets and 3 datasets part of our program) were extracted 
for each source and fit simultaneously. Each dataset was individually reprocessed using CIAO 3.2 and calibration files available in 
CALDB 3.0. 

3. Photometry 

3.1. Optical and near-infrared 



As our objective is to study EROs we proceed by identifying sources in the ISAAC Ks-band image using SEXTRACTOR2.2.2 (Bertin 
& Arnouts|l996Fl and requiring a 3<x threshold above the background in at least 4 pixels. Note that our earlier work on these clus 



ters, including the identification of optical drop-out EROs, was done on an //-band selected sample ( |Richard et al. 2006 1. In contrast 
to Richard et al.|([2006 ) we use AUTOJV1AG instead of aperture photometry, mainly because some of our EROs are quite extended 



and using large enough apertures might induce additional problems due to the close proximity of other sources. As a consequence 



we used the error provided by SEXTRACTOR and not one based on the S/N characteristics in a fixed aperture (see Richard et al. 
2006 ). However, comparing the photometric errors for the EROs already described in Richard et al. (2006 ), we found no significant 



difference between the photometric error based on aperture photometry or AUTO JVIAG. For the SED fitting a minimum photomet 
ric error of 0. 1 mag was generally used. 

For the astrometry we used standard stars from the ESO-USNO-A2.0 catalogue to obtain correct coordinates (J2000). All coordi- 
nates are based on their position in the Ks-band image. 

The photometry in all ISAAC images (SZJHK) was done in the double-image mode of SEXTRACTOR2.2.2, using the Ks-band 
image as reference. 

For objects which were not detected in a specific band, we substituted the detection limit for a 3cr detection threshold as apparent 
magnitude. Using our ERO criterion of (R - Ks) > 5.6 we automatically compiled a catalogue of ERO candidates. The RJHK 
image of each of these candidates was then examined by eye in order to reject spurious detections, e.g. at the edge of the image 
or candidates blended with another source. This procedure resulted in the identification of 6 EROs in Abell 1835 and 9 in AC 114. 
In Richard et al. ( |2006| we found 7 and 1 resolved objects with extremely red optical to near-infrared colours in Abell 1835 and 



AC1 14 respectively, which qualified as ERO. Five of the objects in Abell 1835 cannot be classified as ERO if we use the 3cr detec- 
tion limit for R-band non-detections. Dismissing the optical drop-out criteria applied in that paper and using a different photometry 



1 for details see http://stsdas.stsci.edu/pydrizzle/multidrizzle 

2 This software is freely available from http://terapix.iapJ 
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Fig. 1. Close-up images of EROs found in Abell 1835. Each of the panels is 20 arcsec across, North is up and East is left. Source 
number 1093 is the NIR-counterpart of the sub-mm source SMMJ14009+0252. " The numbers in brackets refer to Richard et al. 
( (2006l > and |SchaereretaIJ l |2007) . 



(see above) we now find 4 additional EROs in Abell 1835 and 8 in AC1 14. The postage stamps of our objects are shown in Figs. [T] 
and|2] 

As result of the higher spatial resolution of ACS, A1835-#319 appears resolved into two sources. It remains to be seen, whether 
or not we are looking at physically distinct sources, or objects which are gravitational bound. For the adopted cosmology the 
projected distance of 0.213 arcsec of the two components of #319 corresponds to 1.7 kpc, assuming a source redshift of 2.5 as 
indicated from their SED analysis (Schaere r et al.|2007] >. These estimates are based on the angular distance without correction for 
lensing and hence state an upper limit for the distance between the two components. 




Fig. 2. Same as Figlllfor the EROs in AC114. Source 512 lies outside the recorded I -band image. 
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Table 4. IRAC and MIPS fluxes for EROS found in Abell 1835 and AC 114. * The source positions at 5.8 and 8 jum are slightly 
displaced from those at shorter wavelengths. Therefore, there is a possibility that the emission at > 4.5/mi is not related to the ERO 
but to the fainter sources to the east. Objects for which we encountered blending problems are not included in Fig.[5]and[6] 



ID 




3.6yum 


4.5fim 


5.8/um 


8.0yum 


24yum 


Abell 1835 




yuly 


/Jy 


fjiy 


yuly 


yuly 


305 


(1) 


2.9± 0.2 


4.0± 0.2 


< 3.6 


<4.5 


< 30. 


319 




19.7+ 0.4 


22.5± 0.4 


26.5± 1.7 


18. 1± 1.6 


272+9* 


347 




15.7± 0.3 


13.9± 0.3 


< 3.6 


8.0± 1.5 


<30. 


532 




73.6± 0.3 


75.3± 0.3 


51.3± 1.5 


37.1+ 1.5 


< 30. 


676 






blended 






1093 


(2) 


14.4± 0.2 


23.0+ 0.3 


37.6± 1.5 


50.9± 1.6 


320+11 


311 


(17) 


2.4± 0.2 


1.6+ 0.2 


< 3.6 




< 30. 


314 


(11) 




blended 






454 


(10) 




blended 






493 


(3) 




blended 






504 


(4) 


1.9± 0.2 


1.0± 0.3 


< 3.6 




< 30. 



AC 114 

572 5.2+0.2 5.3+0.2 7.1 + 1.2 <45 <30T 

572 56.0+0.3 57.4+0.3 51.0+1.3 43.8+1.4 124.9+7.9 

632 210.8+1.0 194.6+0.9 131.2+1.7 92.1+2.2 <30. 

680 27.0+0.2 29.9+0.3 31.8+1.1 16.9+1.4 <30. 

707 46.2+0.3 43.8+0.3 25.1+1.2 23.2+1.4 <30. 

862 39.0+0.2 35.2+0.3 28.7+1.1 22.5+1.4 58.0+10.4 

1006 blended 36.0+1.2 27.5+1.4 51.0+7.9 

1087 48.2+0.2 55.3+0.3 44.0+1.2 36.5+1.5 183.7+7.4 



1167 (1) 67.3+0.5 64.4+0.5 50.3+1.5 44.9+2.2 189.0+8.9 



3.2. IRAC and MIPS photometry 

The IRAC photometry used a circular aperture with a radius of 2'.'4 with a sky background annulus of 2'.'4—7'.'2 in radius. The 
point-source aperture corrections were applied, which were 1.213, 1.234, 1.379, and 1.584 at 3.6, 4.5, 5.8, and 8.0 fim, respectively, 
based on the IRAC Data Handbook. 

The MIPS 24 /vm photometry used a circular aperture with a radius of 6" with a sky background annulus of 6"-13" in radius. 
The corresponding point-source aperture correction was 1.698 based on the MIPS instrument Web site. For several EROs we 
encountered severe blending problems (see Tabj4|i and hence these were not included in our discussion. 



3.3. Chandra photometry 

None of the ERO sources were detected in the fields of AC114 and A1835. In order to place upper limits on the possible X-ray 
emission from the EROs, we have compared the extracted spectra discussed in Section 2.3 at each ERO source position with a 
number of spectral models for the assumed underlying spectral energy distribution. These flux distributions were fit to a simple 
power law model including foreground Galactic absorption. Such a model would be expected if the intrinsic X-ray spectra of the 
ERO was dominated by AGN emission. The absorbing column was fixed to the Galactic value for each cluster and held fixed during 
the fitting procedure. Due to the low number of counts associated with a given source, the spectral index of the power-law model 
was similarly held fixed during the fitting. Values of 1.0, 1.4, and 2.0 were considered for the photon spectral index and different 
energy ranges were considered for comparison with other data from the literature. The resulting flux limits for each ERO source (at 
a 3-sigma level) are listed in Table [5] 



4. Empirical properties of our EROs and comparisons with other samples 

As mentioned above our search for EROs with R - Ks > 5.6 in Abell 1835 and AC114 has yielded 15 (16) objects in total, 
depending if or not sources appearing as double in the ACS images are counted. Two of the additional EROs (#347 & #532) 
detected in Abell 1835 were not included in Ric hard et al.| ( |2006| ) and |Sch aerer et al. (2007 ) due to the additional selection criterion, 
optical non-detection, imposed in these papers. Two more objects (#319 & #676) were previously excluded during the by eye 
examination (as described in section[3]). 

Six of the additional EROs in AC1 14 can be detected in R using the initial detection threshold of lcr and hence were not included 
in Richard et al. ( |2006| l or |Schaerer et al.| ( |2007| l. The remaining two new sources were previously excluded during the visual 
examination, either due to a close by bright object (#1006) or because its position at the edge of the image (#1087). 
We now discuss the properties of the EROs and compare them to other samples and to related objects. At this point we note that we 
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Table 5. X-ray flux limits for the EROs in Abell 1835 and AC 114 computed for two different energy bands and different values of 
the photon energy index F. All upper limits are in units of ergs s^'crrT 2 and refer to a 3cr detection. 



ID 


(0.5-7.0 keV) 


(2.0-10.0 keV) 


Abell 1835 


r=i.o 


T=2.0 


T=1.4 


305 


3.14e-16 


2.21e-16 


2.64e-16 


319 


5.44e-16 


3.85e-16 


4.60e-16 


347 


9.98e-16 


7.08e-16 


8.43e-16 


532 


3.81e-16 


2.70e-16 


3.21e-16 


676 


1.71e-15 


1.22e-15 


1.45e-15 


1093 


5.47e-16 


3.87e-16 


4.63e-16 


311 


5.42e-16 


3.85e-16 


4.60e-16 


314 


3.72e-16 


2.64e-16 


3.16e-16 


454 


3.98e-16 


2.81e-16 


3.36e-16 


493 


4.72e-16 


3.34e-16 


3.99e-16 


504 


5.11e-16 


3.61e-16 


4.30e-16 



AC114 


512 


8.79e-16 


4.52e-16 


6.55e-16 


572 


2.37e-15 


5.75e-16 


1.26e-15 


632 


3.25e-15 


9.74e-16 


2.06e-15 


680 


4.32e-15 


1.70e-15 


3.05e-15 


707 


1.09e-15 


5.58e-16 


8.13e-16 


862 


3.05e-15 


1.55e-15 


2.26e-15 


1006 


1.25e-15 


6.31e-16 


9.25e-16 


1087 


9.72e-16 


4.91e-16 


7.19e-16 


1167 


8.63e-16 


4.38e-16 


6.40e-16 



do not differentiate between the various K-band filters. 



4.1. EROs with very red R - Ks colours 

Fig.[3]and[4]show some of the photometric properties of our EROs in comparison to other samples. These figures show that there are 
some sources, both in Abell 1835 and AC1 14, which have very red colours both with (R - Ks) k 7 and (/ - Ks) ^ 6 (adopting 3cr 
limits for the non-detected bands). Such relatively bright (Ks < 20.3) and red sources were not found in other surveys, e.g. HUDF 
( jYanetal.|2004) l, MUNICS ( |Longhetti et al.|2005|> and GOO DS-MUSIC ( jGrazian et al.|2006|[2007l > although their depth is sufficient 



to detect such red objects. However, Sawicki et al.| ([2005) report 5 EROs with R - K s >7.0 at similar magnitudes as our objects. 
The origin of these apparent differences is not clear. However, it has to be recognised that all samples except GOODS-MUSIC are 
relatively small and lack statistical significance. 



4.2. ERO classification (starburst vs. old population) 

Although colours alone cannot provide the same strong constrains on their nature and photometric redshift as SEDs or spectra, we 
used the colour based classification scheme introduced by Pozzetti & M annucci] ( |2000| l for a first classification of our sample. This 
photometric method uses the (R - K) vs. (J - K) colour plane to separate between galaxies with old stellar populations and dusty 
starbursts, assuming a redshift range of 1 ^ z ^ 2. The corresponding colour-colour diagram of our EROs is shown in Fig. [3] For 
comparison we have also included EROs found by the MUNICS survey (Longh etti et al.| [2005 ), IEROs in the HUDF ( jYan et ah] 
2004) and a sample of EROs found by |Takata efaL] ( p003) in the field of the submillimeter source SMM J04542-0301 (cluster 



MS045 1.6-0305). 

However, taking other information into account this simple classification scheme does not always yield consistent results. For 
example, several objects classified as "elliptical" on the basis of Fig. [3] are detected at 24 //m - incompatible with an old and dust 
free population. Furthermore the SED analysis (cf. Sect. [5]l of these objects and several other objects in the "ellipticals" region 
shows that they are more likely dusty bursting objects. Nevertheless, one has to keep in mind that many of our objects have 



colours close to the separation line (as calculated by |Pozzetti & Mannucci 2000 1 and that the gap between both populations is 
approximately 0.3 magnitudes wide (shaded area in Fig|3|. In this respect, many of our less extreme EROs could also be classified 
as starburst. Also, Pozzetti & Mannucci include exponential declining SFH up to decay times of r ~ 0.3 Gyr in their models of 
evolved populations, while we only consider instantaneous burst scenarios for evolved populations. 
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on a 3cr non- detection (see Rich ard et al.|2 006 ; Scha erer et al.|2007] l. Black diamonds show EROs found by Takata et al. ( 2003| > in 
the vicinity of SCUBA source SMM J0452-0301. Black squ ares mark IEROs in the HUDF ( |Yan et al.|2 004 ) and blue triangles show 
EROs found in the MUNICS survey (Long hetti et al.|2005 1.The dashed line shows the separation between old passive galaxies and 
dusty starburst, according to Pozzetti & Mannucci (2000). The shaded area represents the gap between the ERO populations («0.3 
mag). Arrows indicate the R - Ks or J - Ks colour to be a lower limit due to the non-detection in R or J respectively (see Tabs. [2] 
and|). 
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Fig. 4. I - K vs K colour-magnitude diagram for our sample of EROs (red), IEROs in the HUDF (bla ck, 
EROs found in GOO DS-MUSIC (green, |Grazian et al.|2007) . Blue stars identify the objects described in |Richard et al.| ( [2006| > and 
Schaerer et al. (2007). Arrows mark upper limits, open circles denote EROs without R-band detection. 



2004) and 



4.3. Redshift estimates for EROs with IRAC and/or MIPS detection 



According to Wilson et al. (2004) an ERO selection equivalent to R - K > 5 (less red than our colour threshold) using the IRAC 
3.6/im band would require an R - [3.6] colour redder than 6.6 (Vega) or 4.0 (AB magnitude). All our IRAC detected EROs fulfil 
this criterion (generally R - [3.6]^ fl > 7.4). The combination of Ks-band and IRAC bands also allow a rough estimate of redshift, 
based on the shift of the 1.6 fim bump. This spectral feature can be found in the spectra of all galaxies with the exception of AGN 
dominated SEDs, and hence it can be used to estimate the photometric redshift ( |Sawicki|2002| ). In practice using colour criteria of 
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(K-[3.6]) Vega > 0.9and([3.6]-[4.5])v egfl < 0.47 limits the photometric redshift interval to 0.6 < z < 1.3, while (K-[3 .6])v ega > 0.9 



and ([3.6] - [4.5])y ega > 0.47 should select galaxies with redshifts above 1.3, according to | Wilson et al. ((2004). The corresponding 
colour-colour plot snowing our objects and comparison samples is given in Fig(5] All objects except #305 satisfy this K - [3.6] 
criterion. From their red ([3.6] - [4.5]) we expect that #1093 is clearly above z > 1.3. In fact, detailed SED modelling yields 



photometric redshifts estimates of z p h t ~ 2.0 and 2.8-3. respectively (see Section|5|and Schaerer et al. (|2007[l) in agr eement with 
this simple criterion. However, for the remaining objects with ([3.6] - [4.5]) close to the limit proposed by Wilso n et al.| l |2004) , there 
is only partial agreement for the separation of sources above or below z — 1.3 using the two methods, as can be seen by comparison 
with Table|7] More details on the photometric redshift determination of the EROs are given in Sect|5]and Schaer er et al.| ( [2007] >. 
From a study of 24-jum selected objects Magliocchetti et al. ([2007) suggest that objects with extreme 24-/zm to TJ-band ratios of 
logF24/F/{ 5 -3 are likely z ~ 1.6-3. All 6 MIPS detected objects, except # 572 fall in this category. For #1093 both the colour 
criteria discussed above and SED modelling agree with this classification. For the 3 remaining objects (# 1006, # 1087, # 1167) 
a more complete SED fitting yields, however, photometric redshifts between ~ 0.9 and 1.5. Such a simple criterion may thus 
overestimate the redshift of sources with extreme IR/optical flux ratios. However, since our objects have quite faint MIPS fluxes, 
below the levels of ~ 0.3 mJy discussed by Magliocchetti et al. (2007 1 and Houck et al. (2005) their criterion may be correct for 
more luminous sources. 



4.4. Comparison of EROs with other galaxy populations 

A significant overlap between different galaxy populations selected e.g. according to ERO, DRG, and other criteria is known to 
exist. In our case, e.g. 6 of our 15 EROs satisfy also the DRG selection criterion J - K > 2.3 as indicated in Tables [2] and [3] while 
1 1 objects show a large 3.6//m/zg5o flux ratio used by Yan et al. (2004) to classify IRAC selected objects (IEROs). 



3.0 



2.5 



0.5 



Z<1.J 



o z>1.3 



■» A o a 
■1 % o * 




O305 



0.0 



0.5 
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Schaerer et aL|2007)(blue 



Fig. 5. A'- [3.6]vi.[3.6] - [4.5] colour-colour diagram for extremely red galaxies. The red symbols represent our work. Filled: EROs 
with R-band detection, open: no R-band detection. We also includes the ERos from (Richa rd et al.|2006 
stars). The blue circles denote EROs selected on the basis of their red (R - [3.6])v ega > 6.6 or (K - [3.6 
et al.|2004) . Green and orange symbols show distant red galaxies (DRG) with J 
l 2005| ) respectively. Black squares indicate IEROs by Yan et al. ( 2004| ). 



)vena >l-6) co lour (|Wilson 
K >2.3 from Papovich| (|2006|) and|Labbe et al 



Considering e.g. the near-IR and IRAC colours shown in Fig. 5 , we do not find a distinct difference to other populations: EROs 
( |Wilson et al.|2004| ), DRGs ( |Papovich|20"06] ) and IEROs ( |Yan et al.|2004) . However, from this and from Fig. git is clear that all 
sources without iJ-band detection lie at the outer regions in these plots, indicating somewhat more extreme colours than DRGs, 
which are however shared by some IEROs. Overall most of them correspond to objects with very strong extinction as obtained from 
the SED analysis in Section [5] 

In Fig. [6] we show the I - K versus K - 4.5/um colour, which Labbe et al. 2005 1 u se to separate distant red galaxies (DRGs) 
from z ~2.5 Lyman break galaxies. All samples, DRGs (Labbe et al. (2005 ); Papovich ( 2006[ ), orange and green symbols), IEROS 



( |Yan et aL] ( |2004| ), black symbols) and our EROs (red) occupy the same colour space, with EROs without 7?-band detection again 
occupying the outer regions. We also include the ERO sample of Wilson et al. (2004). Three of the DRGs by Labbe et al. (2005 ) 
are thought to have old stellar populations (orange stars), due to their very red I — K colour. The 5 EROs in our sample, which 
have similar or even redder I — K colours and comparable K - [4.5] colours, were classified as "old evolved" galaxies due to their 
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Fig. 6. Observed / - Kvs.K - 4-.5/im colour-colour diagram. Equal symbols to Fig. [5] Arrows indicate I-band non-detections, open 
circles non-detections in R. 



blue J - K colour. However, from SED fitting (see below) we find that the majority of them are best fit with GRASIL spectral 
templates of very dusty star forming galaxies, which is also supported by their detection at 24 /zm. Such extreme templates were 
not considered by Labbe et al. (2005 ). This shows that not all objects with such extreme I — K colours are "old and dead" galaxies, 
as suggested by Labbe et al. (2005)17 In |Schaerer et al. ( |2007] > we have also shown that the bulk of the IEROs of |Yan et~aL (2004) 
are more likely dusty starbursts than old composite stellar populations. Detailed SED analysis including deep mid-IR observations 
may thus be needed to determine accurately the fraction of "old and dead" galaxies among red distant galaxies, as also pointed out 
by |Kriek et ai] ( [2006| . 



4.5. Starburst vs. AGN classification from X-rays 

None of our EROs is detected in X-rays above the background of the galaxy cluster emission. However, thanks to the depth of the 
observations and the location of the objects away from the maximum cluster emission, the upper limits provide some interesting 
information on the nature of the EROs. 

In Fig. [7] we plot the optical (R-band) flux or upper limits of our objects versus their X-ray flux limits. Also shown are other 
ERO samples from the literature (Mi gnoli et al.|2004 Brusa et al. 2005 Severgnini et al.|2006 1 and curves of constant X-ray to 
optical flux ratios; the range between Fx/Fr = 0.1 and 10 is typical for unobscured type 1 AGN. The majority of our objects 
detected in R are close to the border or outside of this region, indicating that they are likely not unobscured AGN. However, 
deeper X-ray observations might position our sources in the part of the plot occupied by those objects. For the other objects this 
comparison relying on the /?-band is not a strong constraint on their nature. More interesting is the comparison of the MIPS 24 
yum flux with the X-rays, which is very useful to compare AGN or starburst dominated objects, as shown e.g. by |Alonso-Herrero| 
et al.| (2004 1. These data are shown in Fig. [8] together with th e regions of typical X-ray/mid-IR fluxes of local hard X-ray selected 
AGNs and local starbursts taken from Alonso-Herrero et al. ([2004). Clearly the majority of 24 fim detected sources , i.e. 5 out of 
7, have X-ray limits excluding AGNs and compatible with expectations from local starbursts. For the remaining objects we cannot 
conclude firmly on their nature. However, no signature of an AGN is detected. In fact if some of these objects turned out to be 
AGN they would correspond to very faint AGN, given their reasonably well established redshift. For our MIPS detected objects the 
relatively small 8/24 /mi flux ratio is also compatible with starburst dominated objects (cf. Magl iocchetti et aL]2 007 ). We conclude 
that the bulk of our EROs are more likely starburst than AGN dominated at near- to mid-IR wavelengths. 



4.6. Magnification 

Our only selection criteria is based on R - K colour which is not influenced by the magnification effects. However, in order to 
compare the number counts with other surveys, either in the field of lensing clusters (e.g. Smith et al. 2002| l or large surveys 
( |Simpson et al. 2006 ; D addi et al.|2000||Smail et al.|20 02) one has to correct for the magnification of the source flux and the dilution 
of the source plane. Magnification maps were derived following the procedure described in Richard et al. ( 2006) using the mass 
models of Abell 1835 (similar to Smith et al. 2005) and AC1 14 ( |Natarajan et al.|1998[|Campusano et al.|2001) . Given the position 
of each object (in terms of RA and DEC), we then determined the magnification factors. Table|6jlists the necessary correction if the 



A. Hempel et al.: EROs found behind lensing clusters 



13 



28 



27 



26 



25 



24 



23 



22 
10 



o oo * 



-16 




10 

f (2-10 keV) 



Fig. 7. R magnitude vs. 2-10 keV flux (in units of erg cm _2 s _1 ) for our sample (red symbols, see Tab.|iJ, and for other X-ray emitting 
EROs. As none of our sources (red and blue symbols) was detected with Chandra, their Xray flux represents an upper limit. For a 
better clarity of the plot we have obmitted the arrows in the horizontal direction. Triangles (black): Mignol fet al.| ([20Q4), diamonds 



(black): |Brusa et al.| (|2005), small circles (black ) : [Alexander etaL](p003), stars: |Roche et al.| ( |2003| l, large"doT XB S J02 16-0435 



Severgnini et al. 2006; Delia Ceca et al. 2004|l, open squares: Mainieri et al. ( |2002| l. The additional objects from Richard et al 



2006 ) and Schaere r et al.|(|2007|l are ind icated with blue stars. The two dashed lines define the region where unobscured type 1 
AGNs typically lie (see |Fiore et al.|200"3"t|Maccacaro et al.| 1988] >. 



source plane lies at redshift 0.5, 1.0, 1.5, 2.0, 3.0 and 7.0 respectively. As the field of Abell 1835 is positioned away from the cluster 
centre, most of the sources experience a relative small correction. These magnification factors ft also have to be taken into account 
to compute absolute quantities, such as the stellar, SFR, etc. derived in Section [5] 



Table 6. Magnification factors ji from the lensing models of Abell 1835 and AC 114 predicted for various source redshifts z s - The 
values of fj. are dimensionless magnification factors, and not in magnitudes. The magnification factors derived for z s = 1-5 were 
used for the calculation of corrected number counts. 



ID 


z s = 0.5 


z, = 1.0 


z s = 1.5 


z s = 2.0 


z s = 3.0 


z s = 7 


Abell 1835 














305 (1) 


1.13 


1.21 


1.24 


1.26 


1.28 


1.30 


319 


1.15 


1.26 


1.29 


1.31 


1.33 


1.36 


347 


1.17 


1.28 


1.32 


1.34 


1.36 


1.39 


532 


1.33 


1.60 


1.71 


1.77 


1.83 


1.91 


676 


1.34 


1.62 


1.75 


1.81 


1.89 


1.98 


1093 (2) 


1.43 


1.82 


1.99 


2.09 


2.20 


2.34 


311 (17) 


1.17 


1.28 


1.32 


1.34 


1.36 


1.39 


314 (11) 


1.15 


1.26 


1.29 


1.31 


1.33 


1.36 


454 (10) 


1.22 


1.38 


1.44 


1.47 


1.50 


1.54 


493 (3) 


1.27 


1.48 


1.57 


1.61 


1.66 


1.72 


504 (4) 


1.29 


1.51 


1.60 


1.64 


1.70 


1.76 


AC 114 














512 


1.60 


2.74 


3.47 


3.97 


4.61 


5.63 


572 


1.32 


1.71 


1.88 


1.97 


2.06 


2.19 


632 


1.70 


3.03 


3.84 


4.39 


5.06 


6.07 


680 


1.65 


3.03 


4.01 


4.74 


5.75 


7.56 


707 


1.60 


2.71 


3.40 


3.86 


4.45 


5.36 


862 


1.31 


1.69 


1.85 


1.94 


2.04 


2.16 


1006 


1.50 


2.26 


2.65 


2.88 


3.14 


3.50 


1087 


1.20 


1.42 


1.51 


1.56 


1.61 


1.66 


1167 (1) 


1.24 


1.53 


1.64 


1.71 


1.77 


1.86 
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Fig. 8. The 24/jm flux vs. 2-10 keV X-ray flux (red: our work, filled: with /?-band detection, open: no 7?-band detection). Arrows 
assign upper limits, both in X-ray and MIPs flux. Black symbols show X-ray selected sources with 24/im counterparts in the CDF-S 
( Alonso-Herrero et al. 2006 1 . Open black circles show sources not detected in hard X-rays. The Chandra 2-8 keV fluxes have 



been converted to 2-10 keV fluxes assuming a power-law with photon index T=1.4. The blue shaded area is the extrapolation of 
the median hard X-ray-to-mid-IR ratios of local (z<0.12) hard X-ray selected AGNs with mid-IR emission ( |Piccinotti et al.|19 82). 
The purple area is the extrapolation of local starburst galaxies from Ranalli et al. ( |2003 1. The extrapolations were taken from 
Alonso-Herrero et al. ( |2004l [2006) 1 . 



4. 7. Surface density 

After correction for lensing and incompleteness we obtain the surface density of EROs shown in Fig[9] Although we use a slightly 
redder colour threshold than most authors, we compare our cumulative number counts with various other surveys. A comparable 
ERO sample in terms of environment is available from Smith et al. (2002), who studied 10 massive galaxy clusters at z~0.2, 
including Abell 1835. However, the single ERO detected there has a colour of R - K=5.3 and hence does not satisfy our colour 
threshold, while even the brightest of our EROs in Abell 1835 (#532 K=17.97) is not part of this sample. The explanation lies 
both in the smaller field size of the UKIRT camera (1.5 arcmin) and the centering of the field on the central cluster galaxy. As a 
result EROJ 140057+0252.4 (see Smith et al. 2002 ) lies at the very edge of our image, which we excluded due to the low signal-to- 
noise. The bright incompleteness limits (80% at 20.6 mag ) in Smith et al. (2002) might be responsible for their non-detection of our 
fainter objects. The largest ERO survey, using the UKIDSS Ultra Deep survey, has been released recently by Simpson et al. (2006 ). 
However, for now a comparison is only possible at brighter magnitudes, K<20. 1 . 

Fig. [9] show s the cumulative surface density of R - Ks >5.6 EROs in comparison with the samples of Smith et al. ( 2002[ ) and 
Simpson et al^pPOS) (R - K > 5.3 and R - K > 6). 

Up to K=2Q.5 mag , our number counts are slightly lower than those found by ( [Smith et al.|2002 Smail et al.|20 02 ) but agree well 
within the l<x error bars. Our lower number counts are also to be expected due to the redder colour threshold. We estimate the 
cumulative surface density of EROs at Ks <20.5 with (0.97+0.31) arcmin -2 , compared to (1.16+0.17) arcmin -2 and (0.50+0.11) 



arcmin by|Smith et al. (2002 1 for R - K >5.3 and R - K >6.0 (up to the same K-limit). The number counts increase only slightly 



for fainter magnitudes, up to (1.36+0.36) arcmin 2 at Ks <22.0. 

The slope of N(> K)= \0 aK for our sample is »=0.78±0.03 for K<20.0, compared to a=1.05+0.05 from|Daddi et al.|(|2000| and 



a=1.04+0.05 for |Smith et al. ( 2002| ). This slope decreases for fainter magnitudes (20 < K < 22) to a=0. 11+0.01. Such a break in 



the cumulative number counts at K~ 19-20 has been observed in various surveys ( Smith et al.|2002 



the actual values for a are found to be larger by a factor of ~2-3 in ( Smith et al 



|2002[ IDaddiefa 



Daddi et al.|2000l ), although 



2000 ) at fainter magnitudes 



(K >20.5). This flattening could be caused by the absence of evolved ellipticals with fainter magnitudes. 

The preliminary results from UKIDSS UDS EDR ([Simpson et al.|2006] ) (K<20.15) show much higher number counts than any other 
survey (e.g. Daddi et al.|2000 Smail et al.|2002[ |Smith et al.|2002[ )~The authors attribute this result partly to the use of different 
filters and apertures. 



5. Analysis of the observed SED 

5.1. SED fitting method 



To analyse quantitatively the observed SEDs we follow the procedures outlined in 
in |Schaerer et al.| ( |2007| ). We shall only briefly summarise the main points here. 



Schaerer & Pello 



(2005) and described in detail 
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Fig. 9. Cumulative surface density of our sample {R - Ks > 5.6) after correcting for amplification (assuming a source plane at 
z,=1.5) and incompleteness. We also include the ERO sample (R - K > 5.3 and 6.0) of |Smith et al.| ( [2002| found in the fields of 
10 massive galaxy cluster lenses and EROs found by Simpson et al. (2006 1 in the UKIDSS Ultra Deep Survey Early Data Release 
(R - K > 5.3). 



The photometry in all bands except 24 /jm with MIPS has been used. To take uncertainties in absolute flux calibrations between 
different instruments into account we adopt a minimum photometric error of 0.1 mag in most computations. We have used an 
updated version of the Hyperz code from |Bolzonella et al.| ( |2000] > to constrain the redshift, stellar population properties (age, star 
formation history), and extinction of the galaxies studied in this paper. To do so we use in particular a large library of synthetic, 
empirical and semi-empirical spectral templates. The templates are gathered into the 4 following groups: 

1 ) Bruzual & Chariot plus |Cole man et aL ( 1980|l empirical template s galaxies of all Hubble types (hereafter named BC or BCCWW 



2) Starburst SEDs from Schaerer (2002 



group). The theoretical Bruzual models (Bruzual & Chariot 2003), taken here for solar metallicity, include various star formation 
histories representative of different Hubble types. The IMF adopted in these models is the Miller-Scalo IMF from 0.1 to 125 M . 

2003) 1 models at different metallicities extended up to ages of 1 Gyr and considering 



instantaneous bursts or constant star formation (hereafter s04gyr group). These models assume a Salpeter IMF from 1 to 100 M Q . 
3) Empirical or semi-empirical starburst, ULIRG and QSO templates. We use starburst templates from the Calzetti et al. ( 1994) 1 
and |Kinney et al. (1996) atlas and the HST QSO template of Zheng et al. ( 1997 1. To include also more obscured objects we have 
added UV to mill meter band templates of EROs, ULIRGS, starburst and normal g alaxies (HR 10, Ar p 220, M82, NGC 6090, M5 1 , 
M100, NGC 6949) from fits of GRASIL models to multi-wavelength observations ( |Silva et al.| ( |1998| ), named GRASIL group). This 
template group will be used in particular to predict mid-IR to sub-mm fluxes, and hence to estimate total bolometric luminosities, 
after fitting the optical to 8 fim part of the spectrum. The main free parameters we consider are: the spectral template (among a 
group), redshift z , and (additional) extinction (A v ) assuming a Calzetti et al. (2000) law. To increase the diversity of empirical 
or semi-empirical templates and to allow for possible deviations from them, we also allow for an additional reddening. From the 
luminosity distance of the object or, if templates generated by evolutionary synthesis models are used, by scaling the template 
SED to the observed absolute fluxes we obtained the absolute scaling for properties such as stellar mass or the star formation rate 
(SFR). In some cases we also use the bolometric luminosity computed from a GRASIL template to determine the SFR. Finally, the 
absolute quantities must also be corrected for the effects of gravitational lensing. The magnification factors listed in Table [6] are 
used for this purpose. 



5.2. Results 



5.2.1. Abell 1835 

SED fits and the derived properties for the optical drop-out object s #305 (1), #311 (17), #314 (11), #454 (10), #493 (3), #504 
(4), and #1093 (2), corresponding to the objects of Ric hard et al.| ( [2006) ) with IDs in brackets, have already been discussed in 
depth in |Schaerer etaL)([2007) . In cas e of #1093 (2), the known sub-mm source SMMJ 14009+0252, this includes also the SCUBA 
measurements by Ivison et al. (2000). We therefore limit the discussion here to the new objects in this field, i.e. #319, #347, #532, 
and #676. A su mmary of their derive d properties is given in Table [7] For completeness and comparison the derived properties of 
the objects from |Schaerer et al. (20O7J are listed at the bottom of this table. Note the stellar masses have been corrected by (1 + z) _1 
to eliminate an error in the absolute scaling found in Sch aerer et al.| ( |2007[ ). In contrast to the degeneracies found for many of the 
objects discussed in |Schaerer et al.| ( |2007) ) the photometric redshifts of the "new" objects show all well defined best fits at low-z 
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redshifts, z ~ 0.9 to 2.5. Object #319, one of the two newly identified / drop-outs, has the highest photometric redshift, zm ~ 2.4— 
2.5, which is well constrained by the "curvature" measured in the IRAC bands due to the stellar peak at 1.6 //m (restframe). Except 
possibly for #319, none of these objects are detected with MIPS at 24 /mi. 

Object #319, is best fitted without extinction and with templates of elliptical galaxies or maximally old simple stellar populations 
(bursts). From SED fitting, this object is thus best characterised as "elliptical" in agreement with its empirical classification (Pig j3]>. 
However, these SEDs are not able to explain the 24 //m flux. 

The three remaining objects all show clear indications for dust, younger ages, and short star formation histories ("bursts"), 
although they would be classified as "ellipticals" according to their (R-K) vs (J-K) colours (see Fig . [3]> . The estimated stellar masses 
of all these objects are between M* x jj. ~ 3. x 10 10 and 10 12 M with small magnification factors fi. 




1 2 
log(wavelength) [jUm] 



Fig. 10. SED fits with GRASIL templates for selected AC114 objects. Black: #572, red: #1006, #1087, yellow: #1067 (highest-z), 
blue: #1 167. Sub-mm data (APEX/LABOCA) for this southern cluster are being aquired. 



5.2.2. AC114 

Using the templates from the BCCWW and s04gyr groups all except one object have best-fit redshifts of Zfi t ~ 1.-2.6 (see Table |7J. 
The exception is the /^-dropout #1087, whose SED shows a very rapid and strong decline between the near-IR and optical bands, 
which is better fit with a Lyman-break than with a Balmer break. Formally its best photometric redshift is thus zg t ~ 7., but given the 
brightness of this object at near-IR and Spitzer bands the high-z solution is extremely unlikely. Using the semi-empirical GRASIL 
templates the best fit is found at Zfit ~ 0.9. This solution also naturally explains the observed 24 /mi flux of this object, as shown 
in Fig. [10] For #862 the best fit with the BCCWW and s04gyr templates is found at z fit ~ 1.1, where GRASIL templates yield a 
somewhat lower redshift of ~ 0.6 retained in the Table. 

For 4 objects (#512, #632, #680, #707) the SED fitting indicates relatively small amounts of extinction (A v $ 0.6). Interestingly 
these correspond precisely to all the objects which are non-detected at 24 /mi. From their star formation history (all burst-like) 
and their low (but not zero?) extinction these objects resemble ellipticals, in agreement with their empirical classification. For the 
brightest of all our EROs, #632, we obtain different solutions with the BC and S04gyr templates respectively, namely Zfi t ~ 1.3 
(2.5), ages of 4.5 (2.6) Gyr, and little or no extinction, A v ~ 0.6 (0.). The estimated stellar masses of these objects are between 

x yU ~ 6. x 10 10 and 10 12 M . 

The remaining objects, also all detected at 24 /mi, show all clear evidence for substantial extinction. In Table [7] we list their 
properties derived from fits with the GRASIL templates. The best fits are found with the M51 or M82 templates requiring, however, 
additional extinction in 5 of 6 cases. The predicted 24 /mi fluxes, not included in the SED modelling, agree well with the observations 
except maybe for 572 whose MIPS flux is somewhat overestimated (see Fig. 10 1. In conclusion, for 572, 862, 1006, 1087, and 1 167 



all indications (SED fitting and MIPS detections) consistently confirm the dusty low-z (~ 0.9 to 2 or maybe 2.5) starburst nature of 
these objects. This shows the limitations of empirical diagrams, which would e.g. clearly classify three of them as ellipticals (cf. 
Figj3]l. As already mentioned above, our conclusion does also not support the explanation of objects with such extreme I-K colours 
(cf. Fig j6]> as "old and dead galaxies" as proposed by Labbe et al. (2005 1. 

The star formation rates estimated for these dusty objects from the bolometric luminosity of the GRASIL model fit are between 
SFR ~ 15 and 120 M yr~', after correction for lensing. Their bolometric luminosities classify them in the range of luminous 
infrared galaxies (LIRG) with Lboi > 10 11 L . 
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Table 7. Derived/estimated properties for EROs in Abell 1835 and AC1 14. Listed are the object ID (col. 1), the photometric redshift 
estimate (col. 2), the extinction (col. 3), the type of the best-fit template (col. 4), the distance modulus corresponding to z p h t (col. 5), 
the absolute Ks-band magnitude non-corrected for lensing (col. 6), the absolute rest-frame Ks-band magnitude non-corrected for 
lensing (col. 7), the estimated stellar mass (from scaling the SED fit or from M rest (Ks) assuming Lk/M = 3.2, col. 8), the estimated 
star formation rate non-corrected for lensing (col. 9), and the age of the stellar population (col. 10). To correct the above mentioned 
absolute quantities for gravitational magnification the appropriate magnification factors listed in Table [6] must be used. The data 
at the bottom of the Table is from |Schaerer et aX] ( |2007| l (Tabs. 5 and 6), except for the stellar mass, which has been corrected by 
(1 + z)~ 1 to eliminate an error in that paper. 
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5.3. Discussion 



As we can see from Table [7] the properties of our EROs span a rather wide range in extinction, stellar age, and stellar mass. The 
properties and their range are quite similar to those determined for DRGs by Forster-Schreiber et al. (2004) and for the IEROs of 
|Yan et aL] ( |2004 1 in |Schaerer et al.| ( |2007[ ). However, in our sample we find some extreme objects in terms of colours, for which 
the SED modelling indicates quite clearly very high extinction (A v ~ 3 and higher), which is not found in the DRG samples of 
Forster-Schreiber et al.|(2004 1 and Papovich ( 2006 1. According to our analysis (see Schaerer et al.|20O7 1 some IEROs of Yan et al. 



( 2004 ) show also such high extinction. 



6. Summary and conclusions 

We have combined new ACS/HST observations, Spitzer IRAC and MIPS guaranteed time observations, and the optical and near-IR 



observations of Richard et al. (20061 of two well-known lensing clusters, Abell 1835 and AC114, to study extremely red galaxies 
(EROs) in these fields. New and archival X-ray observations with ACIS/Chandra have also been obtained for these clusters. 

Using a standard R - K > 5.6 criterion we have found 6 and 9 EROs in Abell 1835 and AC114 respectively. Several (8) of 
these objects are undetected up to the / and/or z&so band, and are hence "optical" drop-out sources. Three of them, already identified 
earlier by Ric hard et al.| ( |2006] l, have been discussed in detail in |Schaerer et al. (2007]). 

We have discussed the empirical properties of these EROs and compared them to other samples in the literature. We have also 
undertaken SED modelling based on a modified version of the Hyperz photometric redshift code and using a large number of spectral 
templates, including also very dusty galaxies. 

The main results can be summarised as follows: 



- Among our EROs we find 3 sources showing quite unusually red colours in R - K and other colours. Few similar objects are 
found e.g. among the samples of IRAC selected IEROs of |Yan et al"] pXH} , EROs of |Sawicki etlilj ( J2005) and the DRGs of 
Labbe et al. (2005). Our source density is compatible with other counts from the literature. 
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After correcting for lens amplification, we estimate a surface density of (0.97+0.31) arcmin~ 2 for EROs with (R - K >5.6) at 
K< 20.5. We observe a significant flattening of the number count at K~20, possibly the result of loosing the contribution of 
bright evolved ellipticals to the overall ERO population. 

According to "empirical" and to Hyperz modelling, the photometric redshifts of most of our sources yield are z ~ 0.7-1.5. Five 
of them are found at higher redshift (z ~ 2.-2.5.). 

According to simple colour-colour diagrams the majority of our objects would be classified as hosting old stellar populations 
("ellipticals"). However, there are clear signs of dusty starbursts for several among them. These objects correspond to the most 
extreme ones in R - K colour. 

We found that some very red DRGs, which would be classified as old and dead galaxies according to other studies (e.g. Labbe 



et al. (2005 1), are rather very dusty starbursts, even (U)LIRGs, as also supported their mid-IR photometry. Estimates of the 
fraction of old and dead galaxies among red galaxies may thus need to be treated with caution. 

As in earlier studies an overlap of different populations is found. Among our 15 EROs six also classify as DRGs (40%). 12 of 
14 EROs (85 %) with available IRAC photometry also fulfil the selection criteria for IRAC selected IEROs of |Yan et aL] ( |2004| l. 
Objects which do not classify as IERO are also not DRGs; the reverse is however not true. SED modelling shows that ~ 40 % 
of the IEROs are luminous or ultra-luminous infrared galaxies ((U)LIRG). 

None of our objects detected at X-rays above the cluster background emission, with upper limits typically of the order of 
~ (3 — 10) x 10~ 16 ergs s~'cirT 2 in the 0.5-7.0 keV band. No indication for AGNs is found, although faint activity cannot be 
excluded for all objects. From mid-IR and X-ray data 5 objects are clearly classified as starbursts. 

Quantitative SED fitting for our objects shows that they cover a fairly wide range in properties, such as extinction, stellar age, 
mass, and SFR. The derived properties are quite similar to those of DRGs and IEROs, except for 5 extreme objects in terms of 



colours, for which a very high extinction (Ay ^ 3) is found. According to our analysis some IEROs of Yan et al.| (2004) show 



also such high extinction (see Schaerer et al.|2007 1. From the SED modelling these 5 EROs are expected to be (U)LIRG, and 



their IR to sub-mm SED is predicted. 
Understanding the links between these different galaxy populations and their evolutionary history remains largely to be done. 
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